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Abstract 

In this paper we address a class of "ultimate" generation experiments for the 
search of CP-violation in neutrino oscillations. Neutrino factories require large mag- 
netized detectors. New generation superbeams or beta-beams need giant detectors. 
The liquid Argon TPC technology has great potentials for both applications. Al- 
though the ICARUS program has demonstrated that this technology is mature, the 
possibility to built a giant liquid argon TPC is viewed by many as a technically im- 
possible and unsafe task. We argue that a giant liquid argon Cerenkov and charge 
Imaging experiment would be an ideal match for a superbeam or a betabeam. Such 
a detector would in addition cover a broad physics program, including the observa- 
tion of atmospheric neutrinos, solar neutrinos, supernova neutrinos, and search for 
proton decays, in addition to the accelerator physics program. We show a potential 
implementation of such a giant LAr detector and argue that it could be technically 
feasible. The possibility to host such a detector in an underground cavern is under 
study. 
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1 Introduction 



How can one experimentally observe the CP-violation in the leptonic sector? 
From the unitary mixing matrix, which can be parameterized as 



-iS 

C12C13 ■S12C13 SisC 



~Si2C23 — Ci2Si3S23e C12C23 — Si2Si3S23e C13S23 

■S12S23 " Ci2Si3C23e*'^ — C12S23 — Si2Si3C23e*'^ C13C23 

(1) 

with Sij = sin 6ij and Qj = cos 6ij, we get the freedom of the complex phase 
(physical only if ^13 7^ !). We know that 

1. the 5-phase can only be observed in an appearance experiment since the 
disappearance is a T-symmetric process; 

2. the effect for antineutrinos should be opposite to neutrinos {S —S); 

3. it should have the expected L/E^, dependence, where E^, is the neutrino 
energy; 

Considering oscillations involving electron and muon flavors, the oscillation 
probability, in the parameterization described above, is: 



^'-13 



sin^ A23S12S13S23 + C12 (sin^ Aissj^sj^ + sin^ A^sj^ (l - (l + sj^) sl^)) 

— -c^3 sin(26'i2)si3 sin(26'23) cos 6 [cos 2A13 — cos 2A23 — 2 cos(26'i2) sin^ A12] 

+ -c^3 sin 5 sin(26'i2)si3 sin(2^^23) [sin 2Ai2 — sin 2A13 + sin 2A23] (2) 

where Aj^ = Am'^i^L/AEi, (in natural units). 

We note that a precise measurement of the oscillation probability 

can yield information of the (5-phase provided that the other oscillation param- 
eters in the expression are known sufficiently accurately. In practice, one can 
introduce a certain number of discriminants ^|2I which are good quantities to 
experimentally search for a non- vanishing phase: 

1. As = P{ve ^ i^;., 5 = +vr/2) - P{v, ^v^,b = 0) 

The discriminant A^ can be used in an experiment where one is comparing 
the measured z/g oscillation probability as a function of the neutrino 
energy E^ with a "Monte-Carlo prediction" of the spectrum in absence 
of 5-phase. 

2. Acp(5) = P(z/e ^ z/^, 6) - P(z7e ^ Z7^, 5) 

The discriminant Acp can be used in an experiment by comparing the 
appearance of negative and positive muons as a function of the neutrino 
energy E,, for z/g and z/g sources. 



1 



3. At{6) = P(z/e ^ - P{u^ ^ u,,6) or At{6) ^ P{u, ^ u^,5) - 

The discriminant can be used in an experiment by comparing the 
appearance of z/^ and z/g (resp. z/^j and z/g) as a function of the neutrino 
energy E^. 

Each of these discriminants have their advantages and disadvantages^]. The 
A(5-method can provide excellent determination of the phase limited only by 
the statistics of accumulated events, in practice, systematic effects will have 
to be carefully kept under control in order to look for a small effect in a 
seen-data versus Monte-Carlo-expected comparison. In addition, the precise 
knowledge of the other oscillation parameters will be important. The Acp is 
quite straight-forward, it suffers, however, from matter effects which "spoil" 
Acp since it involves both neutrinos and antineutrinos, which oscillate differ- 
ently through matter. Hence, the A^p requires a good understanding of the 
effects related to matter. These effects increase with baseline and are maxi- 
mum around the "resonance energy" around c:^ 10 GeV^. The At is the 
theoretically cleanest method, since it does not suffer from the problems of 
AS and Acp- Indeed, a difference in oscillation probabilities between i/g — >■ z/^ 
and z/^ — >■ z/e would be a direct proof for a non-vanishing 5-phase. In addition, 
matter affects both probabilities in a same way, since it involves only neutrinos 
(resp. antineutrinos for Ay). 

2 Guidelines for a CP-violation neutrino os- 
cillation experiment 

A long-baseline neutrino oscillation experiment designed to search for a non- 
vanishing (5-phase is necessarily an "ultimate" (or at least a "phase-II") ex- 
periment since it should 

• be designed to have ample statistics (for a given ^13) to precisely de- 
termine the oscillation probability as a function of the neutrino energy; 
hence its size will depend on the actual value of 613 (and if ^13 > !) 

• have an excellent energy resolution to observe the energy dependence 
of the oscillation probability and help lift degeneracy of the parameters 
governing the neutrino oscillations (see e.g. Ref.pj); 

• be performed at a wide-band neutrino beam to cover enough "oscilla- 
tions" peaks or do "counting" at different neutrino beam energy settings; 

• have the possibility to study neutrinos and antineutrinos ideally sepa- 
rately in order to lift degeneracies (even in the counting mode). 
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Figure 1: The v^. oscillation probability as a function of neutrino energy for different 
values of the b phase for a given set of oscillation parameters and a distance of 730 km. 

Figure ^illustrates qualitatively the fact that a measurement of the oscillation 
probability as a function of energy with good resolution would indeed provide 
direct information on the (5-phase, since this latter introduces a well-defined 
energy dependence of the oscillation probability, which is different from the, 
say, energy dependence introduced by 6^13 alone (when 6 = 0). 

To study these oscillations, one considers three types of neutrino beams 
produced at accelerators: 

• Superbeams: these are conventional pion beams of high intensity, where 
the sign of the focusing can be selected to enhance neutrinos or antineu- 
trinos components 

7r+ ^ + or 71' ^ + (3) 

• BetabeamsjSj: these are radioactive storage rings, where the type of ion 
can be chosen to select electron-neutrinos or electron-antineutrinos 

^A^ ^-^Al3+ue or ^A^ ^+^A(3-Ue (4) 

• Neutrino factories P*: these are muon storage rings, where the charge of 
the muons can be a priori chosen to select neutrinos/antineutrinos of 
given flavors-pairs 

fi~ -> e"z/eZ/^ or -* e^UeUf, (5) 

From the point of view of neutrino oscillation studies, these beams turn out 
to be very complementary. Indeed, the combination of superbeams (vr^) and 
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beta-beams allows to study CP- and T-violation by a direct comparison 
of the results from the various sources: 

|v ^ (P+) V, (Jt+) 

\ (P-) ^ V (JT) 

In a similar way, a neutrino factory (/i''^) can provide all possible combinations 
to access directly CP- and T-violation: 
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However, it goes without saying that these types of beams require very 
different accelerator technologies. Also from a detector point of view, there 
are quite distinct optimizations required for super- and beta-beams on one side 
and neutrino factory on the other, as discussed in the following sections. In 
particular, detectors for a neutrino factory are necessarily magnetized, hence 
this will limit their size to say 10-50 kton, while detectors for super- and beta- 
beams are conceived as giant non-magnetized experiments with sizes of say 
100-1000 kton. 

2.1 Detectors for the super and beta-beams 

Superbeams and betabeams are typically designed with energies in the range 
of the GeV or below. Since the cross-sections are approximately linear with en- 
ergy, detectors must be very massive, typically in the range of 100-1000 kton. 
For the superbeams, a good e/vr^ discrimination is important in order to sup- 
press the neutral current background with a leading 7r°: 

z/A^ 7r° + X (6) 

from the oscillation channel 

CLnd VeN —>■ e -\- X (7) 

In the case of the betabeams, a good yu/vr^ discrimination is important in 
order to suppress the neutral current background with a charged leading vr^: 

i^N ->TT^ + X (8) 
4 



from the oscillation channel 



and u^N /j, + X 



(9) 



For the beta-beam there is a minimum baseline between the source and 
the detector. Indeed, the oscillation is detected through the muon appearance 
channel with has an energy threshold of about 110 MeV. Hence, 



'mm 



— TT 



nOMeV 
1.27Am2 



136 km 



(10) 



for Am^ = 2 X lO'^ eV^. 

Summarizing, the ideal detector at the superbeam or beta-beam should pos- 
sess the following characteristics: 

• Low energy threshold: the detector should be able to detect, recon- 
struct and analyse events with neutrino energies in the GeV and below. 

• Pcirticle identification: the detector should be able to identify and 
measure electrons and muons, and separate them from other hadrons 
(typ. neutral and charged pions). 

• Energy resolution: the incoming neutrino energy Ei, is reconstructed 

as E,y = Ei + Ehadj where E^ is the leading Icpton energy and Ehad is 
the hadronic energy. Hence, detector with better energy resolution will 
reconstruct the parameter of the incoming neutrino better, and therefore 
the oscillation probability^. 

• Isotropy: at low energy final states particles are produced at all angle, 
in particular, the leading lepton can even be produced backward. This 
is increased on a nuclear target by the Fermi motion. The probably 
most efficient way to detect these particles is to build a large neutrino 
detector, isotropic in nature, capable of measuring equally well particles 
at all angles. This is also a necessary condition for studying astrophysical 
sources and look for proton decay. 

We note that a magnetic field is not necessarily mandatory for these experi- 
ments. 

2.2 Detectors for the neutrino factory 

Neutrino factories are typically designed with energies much higher than the 
GeV. Given energies in the range of 10 GeV, detectors are typically in the range 

■^Note that it is often said that Fermi motion spoils energy resolution at low energy. This is not true 
when the final state is measured completely. Fermi motion can introduce a high momentum imbalance 
(up to 200 MeV), however, the energy available is related to the binding energy in the nucleus and is 
hence small compared to the incoming neutrino energy. 
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of 10-50 kton. At high energy, one could independently study the following 
flavor transitions: 



— > z/e ^ appearance (11) 

— z/^ disappearance, same sign muons (12) 

— > z/^ — r~ appearance, high energy nu's (13) 

— > z/e disappearance (14) 

^ fi^ appearance, wrong sign muons (15) 

— >■ i^T ^ appearance, high energy nu's (16) 



plus 6 other charge conjugate processes initiated from decays. The ideal 
neutrino detector should be able to measure these 12 different processes as a 
function of the baseline L and of the neutrino energy E^! 

Of particular interest are the charged current neutrino interactions, since 
they can in principle be used to tag the neutrino flavor and helicity, through 
the detection and identiflcation of the flnal state charged lepton: 

u^N ^ i~ + hadrons P^N ^ £^ + hadrons (17) 

Hence, the ideal detector at the neutrino factory should possess the following 
characteristics: 

• Particle identification: the detector should be able to identify and 
measure the leading charged lepton of the interaction, in order to tag the 
incoming neutrino flavor. 

• Chcirge identiflcation: the sign of the leading lepton charge should 
be measured, since it tags the helicity of the incoming neutrino. The 
detector must necessarily be magnetized. 

• Energy resolution: the incoming neutrino energy is reconstructed 
as Eu = Ei + Efuid, where Ee is the leading lepton energy and Ehad is 
the hadronic energy. Hence, detector with better energy resolution will 
reconstruct the parameter of the incoming neutrino better, and therefore 
the oscillation probability. 

• Isotropy: one might want to perform various similar experiments at dif- 
ferent baselines. The probably most efficient way to achieve this is to 
build a large neutrino detector, isotropic in nature, capable of measuring 
equally well neutrinos from different sources located at different base- 
lines L. Because of the spherical shape of the Earth, sources located at 
different baselines L will reach the detector "from below" at different an- 
gles. Isotropy of reconstruction is also a necessary condition for studying 
astrophysical sources and look for proton decay. 
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3 The liquid Argon Technology 



The technology of the Liquid Argon Time Projection Chamber (LAr TPC), 
first proposed by C. Rubbia in 1977 0, was conceived as a tool for a com- 
pletely uniform imaging with high accuracy of massive volumes. The opera- 
tional principle of the LAr TPC is based on the fact that in highly purified 
LAr ionization tracks can be transported practically undistorted by a uniform 
electric field over macroscopic distances. Imaging is provided by a suitable set 
of electrodes (wires) placed at the end of the drift path continuously sensing 
and recording the signals induced by the drifting electrons. 

Non-destructive read-out of ionization electrons by charge induction allows 
to detect the signal of electrons crossing subsequent wire planes with different 
orientation. This provides several projective views of the same event, hence 
allowing space point reconstruction and precise calorimetric measurement. 

The main technological challenges of this detection technique have been 
recently summarized elsewhere[6 . They mainly consisted in: (1) techniques 
of Argon purification (2) operation of wire chambers in cryogenic liquid and 
without charge amplification (3) extremely low-noise analog electronics (4) 
continuous wave-form recording and digital signal processing. 

The feasibility of the technology has been demonstrated by the extensive 
ICARUS R&D programme, which included ten years of studies on small LAr 
volumes (proof of principle, LAr purification methods, read-out schemes, elec- 
tronics) and five years of studies with several prototypes of increasing mass 
(purification technology, collection of physics events, pattern recognition, long 
duration tests, read-out technology). The largest of these devices had a mass 
of 3 tons of LAr [3 E] and has been operated continuously for more than 
four years, collecting a large sample of cosmic-ray and gamma-source events. 
Furthermore, a smaller device (50 1 of LAr [91) was exposed to the CERN neu- 
trino beam, demonstrating the high recognition capability of the technique for 
neutrino interaction events. 

The realization of the T600 detector (from design to construction) lasted 
about four years and culminated with a full test of the experimental set-up, 
carried out at surface during 2001. This test demonstrated the maturity of 
the project. All technical aspects of the system, namely cryogenics, LAr pu- 
rification, read-out chambers, detection of LAr scintillation light, electronics 
and DAQ had been tested and performed as expected. Statistically significant 
samples of cosmic-ray events (long muon tracks, hadronic interactions (see e.g 
Figure EI), spectacular high-multiplicity muon bundles, electromagnetic and 
hadronic showers, low energy events) were recorded. The analysis of these 
events has allowed the development and fine tuning of the off-line tools for 
the event reconstruction and the extraction of physical quantities. It has also 
demonstrated the performance of the detector in a quantitative way. 
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Figure 2: Hadronic interaction collected during the ICARUS T600 technical run on the 
surface. 

The detector performance can be summarized as: 

• a tracking device with precise event topology reconstruction 

• momentum estimation via multiple scattering 

• measurement of local energy deposition (dE/dx), providing e/vr^ separa- 
tion (sampling typ. 2%Xq), particle identification via range versus dE/dx 
measurement 

• total energy reconstruction of the event from charge integration (the vol- 
ume can be considered as a full-sampling, fully homogenous calorimeter) 
providing excellent accuracy for contained events 

The energy resolutions are: 

• a/E = 11%/ ^E{MeV)® 2% for low energy electrons fmeasured|lUj) 

• a/E ^ 3%/ a/ E{GeV) for electromagnetic showers 

• a/E ^ 30%/ a/ E{GeV) for hadronic showers (pure LAr) 

• a/E ^ 11%/ ^E{GeV) for hadronic showers (TMG doped LAr) 

It is fair to say that the technique has reached a high level of maturity. 

4 A magnetized Liquid Argon TPC 

Liquid argon imaging provides very good tracking with dE/dx measurement, 
and excellent calorimetric performance for contained showers. This allows for 
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a very precise determination of the energy of the particles in an event. This is 
particularly true for electron showers, which energy is very precisely measured. 



B=1 



Figure 3: Magnetized liquid argon TPC: simulation of the 2.5 GeV electron shower in 
liquid argon. The field has a strength B=l T and is directed perpendicular to the sheet- 
plane. 

The possibility to complement these features with those provided by a 
magnetic field has been considered and would open new possibilities: 

• charge discrimination 

• momentum measurement of particles escaping the detector (e.g. high 
energy muons) 

• very precise kinematics, since the measurements are multiple scattering 
dominated (e.g. Ap/p ~ 4% for a track length of L = 12 m and a field 
of B = IT). 

The orientation of the magnetic field is such that the bending direction is in 
the direction of the drift where the best spatial resolution is achieved (e.g. in 
the ICARUS T600 a point resolution of 400 /im was obtained). The magnetic 
field is hence perpendicular to the electric field. The Lorentz angle is expected 
to be very small in liquid (e.g. ~ 30mrad at E = 500 V/cm and B = 0.5T). 
Embedding the volume of argon into a magnetic field would therefore not alter 
the imaging properties of the detector and the measurement of the bending of 
charged hadrons or penetrating muons would allow a precise determination of 
the momentum and a determination of their charge. 

The required magnetic field for charge discrimination for a path x in the 
liquid Argon is given by the bending 

/2 _ 0.3B{T){x{m)Y 
2R 



2p{GeV) 

and the multiple scattering contribution: 

0.02(x(m))3/2 



;i8) 



MS 



p{GeV) 



(19) 
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Figure 4: Magnetized liquid argon TPC: ongoing R&D to test liquid argon imaging in 
a magnetic field, (left) CAD drawing of the chamber; the direction of the fields and a 
bending track is shown (right) actual setup. 



The requirement for a 3cr charge discrimination can be written as: 
2b > SMS, which implies: 

0.2(T) 



B > 



'x(m] 



(20) 



For long penetrating tracks like muons, a field of O.IT allows to discriminate 
the charge for tracks longer than 4 meters. This corresponds for example to 
a muon momentum threshold of 800 MeV/c. Hence, performances are very 
good, even at very low momenta. 

We have recently started studying the effect of the magnetic field on elec- 
trons (see Figure ^ by means of full detector simulation and reconstruction. 
Unlike muons or hadrons, the early showering of electrons makes their charge 
identification difficult. The track length usable for charge discrimination is 
limited to a few radiation lengths after which the showers makes the recogni- 
tion of the parent electron more difficult. In practice, charge discrimination is 
possible for high fields: 



x = lXo^ B > 0.5T 
x = 2Xo^ B > OAT 
x = 3Xo^ B > 0.3T 



(21) 
(22) 
(23) 



From the simulation, we found that the determination of the charge of elec- 
trons of energy in the range between 1 and 5 GeV is feasible with good purity, 
provided the field has a strength in the range of 1 Tesla. Preliminary estimates 
show that these electrons exhibit an average curvature sufficient to have elec- 
tron charge discrimination better than 1% with an efficiency of 20%. Further 
studies are on-going. 
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Figure 5: Artistic view of the LANNDD design [12] . See Ref. for explanation of numbers. 

In parallel, we have initiated an R&D to test liquid argon imaging in a 
magnetic fieldlll^ (See Figure Ej). We have built a small liquid argon TPC 
(width 300 mm, height 150 mm, drift length 150 mm) and placed it in the 
recycled SINDRUM-I magnet*' which allows us to test fields up to 0.5 T. The 
ongoing test program includes (1) checking the basic imaging in B-field (2) 
measuring traversing and stopping muons (3) test charge discrimination (4) 
check Lorentz angle. Results are expected in 2004. 

The design of a magnetized liquid argon TPC of 70 kton has been considered |12j. 
It is based on a direct extrapolation ( "scaling up" ) of the technique developed 
by the ICARUS Collaboration, however, embedded in a very large magnet 
(See Figure El). A magnetized liquid argon detector would offer unequaled 
physics opportunities at a neutrino factory [H [12]. However, the technical fea- 
sibility of the enormous magnet with its gigantic yoke remains an unsolved 
challenge. Further engineering studies are mandatory before this technique 
can be proposed in an experiment of this scale. 

''The magnet was kindly lend to us by PSI, Villigen. 
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Figure 6: Super- and Hyper-Kamiokande (see Ref.pU]). 



5 Water versus Liquid Argon 

Both water Cerenkov and liquid argon detectors will cover a broad physics 
program, including the observation of atmospheric neutrinos, solar neutrinos, 
supernova neutrinos, and search for proton decays, in addition to the acceler- 
ator physics program. 

Giant water detectors in the range of a megaton are perceived as a "straight- 
forward" extrapolation of existing detectors like SuperKamiokande|19j. There 
have been several proposals, as an example the design of the Hyper-Kamiokande [211] 
considered in the context of the JHF program is reported in Figure IHl 

On the other hand, the liquid argon technique is considered as a difficult 
or unsafe, an almost impractical, solution for large detectors. In this section, 
we discuss briefly the physical properties of water and liquid argon. These are 
summarized in Tabled The densities of the two liquids are very similar, with 
a 40% advantage for liquid Argon. The radiation and interaction lengths are 
rather similar, even though liquid Argon has roughly half the radiation length 
than water. The stopping power dE/dx of particles is similar in water and 
liquid Argon. Hence, neutrino interactions or other rare events will develop 
very similarly in liquid argon or water! 

The refractive index in the visible spectrum is 1.33 for water and 1.24 
for liquid Argon. This means that the Cerenkov emission properties of both 
media are very similar. Indeed, the Cerenkov angle 9c is 42° in water and 36° 
in liquid argon. The number of Cerenkov photons produced per track unit 
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Property 


Water 


Liquid Argon 




1 
i 




Radiation length (cm) 




14. U 


Interaction length (cm) 


od.o 


00.0 


dill 1 ax \NiB V / cm j 


1 

i.y 




JTlyCii cLL/Ll Vc liiU-CA. I VlolUiCl 






Cerenkov angle 


42° 


36° 


Cerenkov (fN/dEdx (eV"-^cm"^) 


160 


130 


Muon Cerenkov threshold (MeV/c) 


120 


140 


Scintillation 7/MeK @ E = 


No 


Yes, 5 X 10^ 


Electron mobility 




500 cm^lVs 


Long electron drift path 


no 


possible 



Table 1: Comparison between water and liquid Argon. 



length and per unit photon energy is 

- " ■ ^6* (24) 
dEdx he 

^ 370sin2^c eV^cm"^ (25) 

Accordingly, it is ~ 160 eV^^cm"^ for water and ~ 130 eV"^cm~^ for liq- 
uid Argon. Hence, both liquids have similar Cerenkov imaging capabilities. 
Cerenkov light from penetrating muon tracks has been successfully detected 
in a liquid Argon TPC|15j. 

An advantage of noble liquid gases compared to water is their very high 
scintillation light (luminescence) yield, comparable to that of Nal crystals. 
The average energy needed to produce a photon is defined as W^. Its measured 
value (at zero electric field) is 19.5 eV for liquid Argon, 52 eV for hquid 
Krypton and 38 eV for liquid Xenon. This scintillation is produced [IHl HZj 
by the de-excitation of a dimer molecule [Ar — Ar^)* by the emission of a 
single UV photon (see Figure [7j). In the case of liquid Argon, scintillation 
appears essentially as a monochromatic line at a wavelength A = 128 nm. The 
produced light via this mechanism is not energetic enough to further ionize, 
hence, the medium is transparent to its scintillation. In addition, it should be 
stressed that in the liquid argon technique the medium is purified to very high 
levels (e.g. less than 1 in 10^° oxygen-equivalent impurities) in order to ensure 
drifts of electrons over long distances without impurity attachment losses, 
hence it has excellent light transmission capabilities. The propagation of light 
is dominated by Rayleigh scattering. The scattering length is computed to 
be about 90 cmjTH] for the scintillation line at A ~ 128 nm. Owing to the 
A^ dependence of the Rayleigh scattering process, we expect no effect for 
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Figure 7: Processes induced by charged particles in liquid noble gases (From Ref . pTTI ) ■ 
For particle with velocity (3 > 1/n there is also Cerenkov light emission (not included in 
the chart). 



visible light over distances of lOO's of meters. Visible photons therefore travel 
the medium in straight lines essentially unperturbed. It is also possible to 
increase the apparent scattering length of the scintillation light by doping the 
liquid argon with xenon (typ. 5% mixture) which shifts the light to longer 
wavelengths. 

6 Combining Cerenkov, scintillation light and 
charge readout in a liquid Argon TPC 

Since water and liquid Argon have very similar Cerenkov light emission proper- 
ties and also similar physical properties in terms of radiation lengths, interac- 
tion lengths, etc.. the events in water and in liquid Argon look very much the 
same and the techniques developed in Kamiokande and Sup erkamiokande 1 1 9 j 
for the reconstruction and analysis of events can be readily "transposed" to 
the liquid argon case. Hence, the performance of a liquid argon detector with 
Cerenkov light readout is at least equivalent to that of a giant water Cerenkov 
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detector from the point of view of event detection, reconstruction and analysis. 

Of course, the overall performance of the liquid argon TPC profits greatly 
from its tracking imaging properties not available in Cerenkov imaging. The 
non-destructive, multiple plane readout allows to reconstruct images in space. 
With the imaging quality of a bubble-chamber all particles can be fully recon- 
structed in space, with extremely low thresholds, and the excellent calorimetry 
allows to reconstruct energies very precisely. Tracking and calorimetry provide 
momenta, particle identification, clean e/vr" separation, etc. Figure |H1 shows 
for example the reconstruction of a neutral pion via Cerenkov and charge 
imaging. Multiple showers are difficult identify in water as their Cerenkov 
rings tend to easily overlap. 




Figure 8: Neutral pion reconstruction via (a) Cerenkov (left) (b) charge imaging TPC 
(right) . 

The most striking difference between water Cerenkov and liquid argon 
imaging is illustrated in Table |2l where the Cerenkov emission momentum 
threshold and the corresponding range in liquid argon are listed for various 
particles. For example, a proton becomes visible in water Cerenkov detector 
when its momentum is greater than 1070 MeV. At this momentum, a proton 
has a range of about 80 cm in liquid argon! With a typical wire pitch of 
3 mm, particles like kaon, protons, etc. can be well detected down to very low 
momenta. 

How can one then profit from the readout of the Cerenkov light in addition 
to the imaging? The combination of the information from the tracking, energy 
(e.g. dE/dx and kinetic energy) with the Cerenkov light provides improved 
particle identification. In particular, one can in this way separate pions from 
muons, a very important tool in the context beta-beams, as illustrated in 
Figure IHl 
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Particle 


Cerenkov thr. in H2O (MeV/c) 


Range in LAr (cm) 


e 


0.6 


0.07 




120 


12 


71 


159 


16 


K 


568 


59 


P 


1070 


80 



Table 2: Momentum threshold for Cerenkov light emission and corresponding range 
liquid Argon for various particles. 



Ui^ combinaiion 
of charge imagiug 
{J(dH/dx)dx= Tj^J 
Sc Cerenkov light 
(P) LO separate 
muoiis from pions 



Figure 9: Illustration of signal and background given by pion/muon confusion in bet 
beams, when searching for z/g oscillations. 

The number of photoelectrons between the wavelengths of 160 and 600 nm 
(the typical acceptance of a PMT) emitted by muons and pions along their 
trajectory before they range out is shown in Figure [TUT left) as a function of 
the kinetic energy of the particle. We assumed a 20% coverage and a 20% 
quantum efficiency and that photoelectron counting is possible. In the kinetic 
region of interest, the number of photoelectrons varies from 1000 to 10000. 
Pions produce as expected slightly less photons. The statistical error is so 
precise that it should allow separation between the two hypotheses. For the 
pion, we assumed that no hadronic interaction takes place along its range. In 
reality, the observation of an hadronic interaction via the imaging can be used 
to discriminate pions versus muons. The charged pion survival probability 
for 90% muon acceptance efficiency is plotted in Figure ITUT right). In the 
range between 100 MeV (threshold) and 850 MeV, the method is effective at 
separating muons from pions. 
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Figure 10: (Left) Number of photoelectrons between 160 and 600 nm emitted by muons 
(line) and pions (dashed). (Right) charged pion survival probability for 90% muon ac- 
ceptance efficiency. Both quantities are plotted as a function of the kinetic energy of the 
particle. We assumed a 20% coverage and a 20% quantum efficiency. For the pion, we as- 
sumed that no hadronic interaction takes place along its range. In reality, the observation 
of an interaction via the imaging suppresses further pions versus muons. 

7 A Giant Liquid Argon TPC with charge imag- 
ing, scintillation and Cerenkov light readout 

Although the liquid argon TPC technology has been demonstrated to be ma- 
ture, the possibility to construct a giant liquid argon TPC remains for many 
an impossible technical task. In this section, we describe some issues that to 
our mind show that giant liquid argon detectors might be technically feasible. 

The ICARUS collaboration has proposed an underground modular T3000 
detector for the LNGS laboratory based on the cloning of the T600 detector |21j. 
The T3000 would be composed of T600 + T1200 + T1200. The design is fully 
proven by the successful technical run on the surface. Further modules are 
ready to be built by industry. A 10 kton detector based on this design "could" 
be ordered today, even though this would not be most optimal financially. Fol- 
lowing a successful scaling up strategy, one could envision building bigger su- 
permodules based on the ICARUS-Airliquide technique, by readily increasing 
the dimensions of the planned ICARUS T1200 by a factor 2 in each directions: 
this would yield (2^) x T1200 ~ TIOK. Hence, it seems conceivable to scale 
up the ICARUS dewar to a 10 kton volume. However, to reach a total mass 
of 100 kton would still require a large number of such supermodules (in fact. 
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10 X 10 kton ~ 100 kton, note however that the fiducial volume is less given 
the modularity). The modularity increases the complexity of the system. 

It appears from the above discussion that contrary to a modular approach, 
a single giant volume is the most attractive solution. In fact, it appears that 
the maximum size of the single module is limited by the requirement to locate 
the detector underground in a cavern and not by the possibility to build a 
large cryogenic tanker of the needed size. Is a strong R&D program required 
to extrapolate the liquid Argon TPC to the 100 kton scale? Or can it be 
achieved in say one (or two) step(s)? In the following, we try to address the 
feasibility of a single volume 100 kton liquid argon detector. 

7.1 Overview of the basic design parameters 

An artistic view of the detector is shown in Figure ^2 A summary of param- 
eters are listed in Table El The detector can be mechanically subdivided into 
two parts: (1) the liquid argon tanker and (2) the inner detector instrumen- 
tation. For simplicity, we assume at this stage that the two aspects can be 
decoupled. 




Figure 11: An artistic view of a 100 kton single tanker liquid argon detector. It appears 
that the feasibility of a volume of this size will be limited by the requirement to find a 
geologically stable underground cavern of this size. The electronic crates are located at 
the top of the dewar. 

The basic design parameters can be summarized as follows: 
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1. Single 100 kton "boiling" cryogenic tanker with Argon refrigeration (in 
particular, the cooling is done directly with Argon, e.g. without nitrogen) 

2. Charge imaging + scintillation + Cerenkov light readout for complete 
event information 

3. Charge amphfication to allow for extremely long drifts: the detector is 
running in bi- phase mode. In order to allow for long drift (~ 20 m), we 
consider charge attenuation along drift and compensate this effect with 
charge amplification near anodes located in gas phase. 

4. Absence of magnetic field 



Table 3: Summary parameters of the 100 kton liquid Argon detector 



Dewar 


$ Ri 70 m, height 20 m, passive perlite insu- 
lated, heat input Ri 5 W/m^ 


Argon Storage 


Boiling argon, low pressure (< 100 mb overpres- 
sure) 


Argon total vohimc 


73118 m'^ (height = 19 m) .ratio 

area,/ vohunc~ 15%: 


Argon total mass 


102365 TONS 


Hydrostatic pressure at bottom 


~ 3 atm 


Inner detector dimensions 


Disc $ « 70 m located in gas phase above liquid 
phase 


Electron drift in liquid 


20 m maximum drift, HV= 2MV for 
E=l kV/cm, Vd ~ 2 mm//xs, max drift 
time ~ 10 ms 


Charge readout views 


2 independent perpendicular views, 3 mm pitch, 
in gas phase (electron extraction) with charge 
amplification 


Charge readout channels 


« 100000 


Readout electronics 


100 racks on top of dewar (1000 channels per 
crate) 


Scintillation light readout 


Yes (also for triggering), 1000 immersed 8"PMT 
with WLS (TPB) 


Visible light readout 


Yes (Cerenkov light), 27000 immersed 8"PMTs 
or 20% coverage, single photon counting capa- 
bility 



7.2 The 77'000 liquid argon tanker 

In order to achieve such large volumes of liquid argon, we base our design on 
the large industrial expertise in the storage of liquefied natural gases (LNG, 
T ~ llOi^ at 1 bar). The LNG technology has been developed quite dramati- 
cally in the last decades and was driven by the petrochemical and space rocket 
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industries. The technical problems associated to their design, construction 
and operation have already been addressed and solved by the petro-chemical 
industry. The current state-of-the-art contemplates tankers of 200'000 m^. 
Currently there seem to be in the world about 300 giant cryogenic tankers 
with volumes larger than SO'OOO m^. Large ships transporting volumes up to 
145'000 of LNG often cross the oceans. 

LNG tanks are always of double-wall construction with efficient insulation 
between the walls. Large tanks are low aspect ratio (height to width) and 
cylindrical in design with a domed roof. Storage pressures in these tanks are 
very low. LNG are used when volume is an issue, in particular, for storage. 

105,000 LNG Storagfl Tank 




Figure 12: Engineering design of a 105'OOOm'^ cryogenic tanker developed by Technodyne 
International Limited (see Ref.|22j). 

Of course commercial tankers are located on the surface and hence our 
tanker, although of reasonable size, must face the additional constraint of 
being located underground. We have contacted the Technodyne International 
Limitedj221 UK to initiate a feasibility study in order to understand 

what are the issues related to the operation of a large underground liquid argon 
detector. Technodyne is engineering company specialized in large LNG tankers 
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(See Figure[T21). In the baseline configuration, we are studying the design of the 
standard tanker (similar to those on the surface) to be located underground. 
The tanker will be self-supporting and will not rely on the surface of the 
cavern. Initial considerations seem to indicate that the extrapolation from 
LNG to liquid Argon is rather straight-forward. The geophysics of the cavern 
can be understood and possible movements can be predicted for periods of 
time extending to at least 30 years. The cooling of the cavern due to heat 
losses is also taken into account. 



7.3 The inner detector instrumentation 

A schematic layout of the detector is shown in Figure UHl The detector is 
characterized by the extremely large volume of argon. A cathode located 
near the bottom of the tanker is set at —2MV creating a drift electric field 
of 1 kV/cm over the distance of 20 m. In this field configuration ionization 
electrons are moving upwards while ions are going downward. The electric field 
is delimited on the sides of the tanker by a series of ring electrodes (race-tracks) 
put at the appropriate voltages (voltage divider). The breakdown voltage of 
liquid argon is such that a distance of about 50 cm to the grounded tanker 
volume is electrically safe. For the high voltage we consider two solutions: (1) 
either the HV is brought inside the dewar through an appropriate custom- 
made HV feed-through or (2) a voltage multiplier could be installed inside the 
cold volume. 




Figure 13: Schematic layout of a 100 kton liquid Argon detector. 
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The relevant parameters of the charge readout are summarized in Table |3] 
The tanker contains both liquid and gas argon phases at equilibrium. Since 
purity is a concern for very long drifts of the order of 20 meters, we think that 
the inner detector should be operated in bi-phase mode, namely drift electrons 
produced in the liquid phase are extracted from the liquid into the gas phase 
with the help of an appropriate electric field. Our measurements show that the 
threshold for 100% efficient extraction is about 3 kV/cm. Hence, just below 
and above the hquid two grids define the appropriate liquid extraction field. 



Table 4: Parameters of the charge readout 



Electron drift in liquid 


20 m maximum drift, HV=2 MV for 
E=lkV/cm, Vd — 2mm/ fis, max drift time 

tmax ^ 10 ms 


Charge readout views 


two independent perpendicular views, 3 mm 
pitch 


Maximum charge diffusion 


— 2.8mm {■\/2Dtmax for D = Acrri^ / s) 


Maximum charge attenuation 


^-tmax/r ^ 1/150 for T = 2 ms electron lifetime 


Needed charge amplification 


10^ to 10^ 


Methods for amplification 


Extraction to and amplification in gas phase 


Possible solutions 


Thin wires+pad readout, GEM, LEM, ... 



In order to amplify the extracted charge, one can consider various options 

(1) amplification near thin readout wires (like in the MPWC, see Figure IT^ : 

(2) GEM[2n| or (3) LEM[24J. Generally speaking, amplification is technically 
challenging since one has to operate in pure argon and one has to face prob- 
lems of sparking, instabilities, etc. The addition of a quenching gas is not 
practical since it would pollute the liquid. It would also absorb the scintilla- 
tion light from the liquid argon. Nonetheless, we have experimentally obtained 
encouraging results which show that a gain of 100-1000 is achievable. Since 
the readout is limited to the top of the detector, it is practical to route cables 
out from the top of the dewar where electronics crates can be located around 
the dewar outer edges. 

Amplification operates in proportional mode. After maximum drift of 20 m 
at 1 KV/cm, the electron cloud diffusion reaches approximately 3 mm which 
is the size of the readout pitch. Hence, drifting along longer distances would 
start degrading the quality of the images because of the diffusion of the charge 
on adjacent wires. We thus think that 20 meters corresponds to a maximum 
conceivable drift distance. Drifting over such distances should be possible, 
allowing for some charge attenuation due to attachment to impurities. If we 
assume that the reachable electron lifetime is at least r ~ 2 ms (this is the 
value achieved on the ICARUS T600 detector during the technical run [23] and 
better values up to 10 ms were obtained on smaller prototypes during longer 
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Figure 14: GARFIELD simulation for the amplification near wires in pure argon (Ar 
100%, T=87K, p=l atm). (left) Geometrical setup (right) gain as a function of wire 
diameter. The induced signals are (1) wires and (2) strips provide two perpendicular 
views. 



runs), then we expect an attenuation of a factor ~ 150 over the distance of 
20 m. This loss can be compensated by the proportional gain near the anodes. 

In addition to charge readout, we envision to locate PMT's around the 
tanker. Scintillation and Cerenkov light can be readout essentially indepen- 
dently. One can profit from the ICARUS R&D which has shown that PMTs 
immersed directly in the liquid Argon is possible 6 . One is using commercial 
Electron Tubes 8" PMTs with a photocathode for cold operation and a stan- 
dard glass window. In order to be sensitive to DUV scintillation, the PMT 
are coated with a wavelength shifter (Tetraphenyl-Butadiene) . 

As already mentioned, liquid Argon is a very good scintillator with about 
50000 7/MeV". However, this light is essentially a line at A = 128 nm. 
Cerenkov light from penetrating muon tracks has been successfully detected in 
a liquid Argon TPC^HI- This much weaker Cerenkov light (about 700 ^jMeV 
between 160 nm and 600 nm for an ultrarelativistic muon) can be separately 
identified with PMT's without wavelength shifter coating, since their efficiency 
for the DUV light will be very small. 

Summarizing about 1000 immersed phototubes with WLS would be used 
to identify the (isotropic and bright) scintillation light. While about 27000 
immersed 8" -phototubes without WLS would provide a 20% coverage of the 
surface of the detector. As already mentioned, these latter should have sin- 
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gle photon counting capabilities in order to count the number of Cerenkov 
photons. 

7.4 Cryogenic aspects 

In order to guarantee the safety of a tanker of this size one cannot rely on 
vacuum insulation. On the contrary one can use rather convential heat insu- 
lation like for example thick layers of perlite. We envisioned an effective heat 
input average over the area of the detector at the level of 5 W/m^. This is 
to be compared to vacuum insulations which could reach 0.1 W/m^, however, 
with the risk of a vacuum loss. The tanker envisaged has a very favorable 
ratio of area over volume of about 15%. Hence, the heat input of 60 kW has a 
small effect on the big volume. See Table For comparison, the refrigeration 
system of the CERN LHC has a cooling capacity of 140 kW at 4.5 K and 
CERN has since 2003 a total cryogenic capacity of 162 kW at 4.5 K. 

The insulation, as efficient as it is, will not keep the temperature of LNG 
cold by itself. The liquid is stored as a "boihng cryogen" that is, it is a very 
cold liquid at its boiling point for the pressure it is being stored. The hquid 
will stay at near constant temperature if kept at constant pressure. This 
phenomenon is called "autorefrigeration" . As long as the steam (liquid vapor 
boil off) is allowed to leave the tank, the temperature will remain constant. 
Hence, the safest way to store large quantities of cryogenic liquid is at a small 
overpressure of less than approximately 200 mbar (basically at atmospheric 
pressure) and let it evaporate. While the liquid evaporates, the temperature 
remains constant. Since the process of evaporation is slow, the tanker remains 
in very stable and safe conditions. 



Table 5: Cryogenic parameters: boiling (Heat loss should be conservative for 3 meter 
thick perlite and includes heat input from supports, instrumentation (cables), etc) 



Dewar 


$ « 70 m, height~20 m, passive 3 m thick per- 
lite insulated, assumed effective heat input ~ 5 
W/m2 


Total area 


12100 


Total heat input 


60500 W 


Liquid Argon evaporation rate 


0.27 liters/second or 23000 liters/day 


Fraction of total evaporation rate 


0.03% of total argon volume per day 


Time to totally empty tanker by evap- 


9 years 


oration 
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7.4.1 Initial filling 



Because of the large amount liquid argon needed to fill up the experiment 
(e.g. 300 ton/day to fill in 300 days), we think that liquid argon could be 
produced locally. Other options (transport from various centers) are being 
investigated. If we envision a dedicated cryogenic plant, then it should be 
located on the surface (i.e. not underground) and connected to the detector 
via long vacuum-insulated pipes. See Figure UHl Hot argon gas is compressed 
to high pressure (typ. 200 bar) and transported towards the underground 
complex. In doing so, it can be purified. A heat exchanger is located before a 
standard Joule- Thompson expansion valve which liquifies the argon into the 
cryogenic tanker. The gas phase is extracted from the tanker through a heat 
exchanger and transported outside. Then it can either be ventilated other fed 
into the compressor for a further cooling cycle. This simple process is called 
a Linde-refrigerator. 




Figure 15: Cryogenic setup of the 100 kton liquid Argon detector. 



The advantage of Argon is that it is naturally available everywhere in the 
world, since it composes approximately 1% of air (by volume air is 78,1% 
nitrogen, 20% oxygen, 0.93% argon, 0.033% carbon dioxyde and 0.003% rare 
gases like neon, helium, krypton, hydrogen and xenon). Liquid Argon is hence 
extracted from the standard process of liquefaction from air. Air mixture 
is cooled down and cold gas-mixtures are separated into Oxygen, Nitrogen, 
Argon, ... Liquid Argon is used to fill the experiment. The rest of the gases 
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can be sold'^. An energetically efficient method is to employ the "useless" 
liquids to improve the efficiency of the liquefiers. 

From the energy balance point of view, the Linde-refrigerator and the pro- 
cess of separation of Argon from air can be ideally computed with basic ther- 
modynamics. For a mixture of ideal gases, one has|26j 

|iy,| = i?T^2/,ln(l/i/,) (26) 

3 

where Uj is the mole fraction of the jth mixture. It turns out that the ideal 
power of separation of Argon from air is 354 kJ/kg at T = 300 K. Similarly, 
the ideal liquefaction process requires a work[26] 

\Wi\=m{T{s,-Sf)-{K-hf)) (27) 

where s is the entropy and h the enthalpy of the fluid (initial, final). The ideal 
argon liquefaction is 478 kJ/kg at T = 300/^, p = latm. If we assume an 
overall thermodynamical efficiency of 5%, then we find that the total electrical 
power needed to ffil the 100 kton tanker in two years is 30 MW. Indeed, to ffil 
the tanker in such short time requires 150 tons of liquid argon per day! These 
figures are summarized in Table IHl 



Table 6: Cryogenic parameters: initial filling phase (Initial cooling of tanker not 
included) 



Liquid Argon 1st filling time 


2 years (assumed) 


Liquid Argon 1st filling rate 


1,2 liters/second or 150 tons/day 


Argon gas equivalent 


85000 mVday 


Air volume equivalent (Ar 1%) 


8500000 mVday (205 m)Vday 


Ideal power of separation of Argon 


600 kW (assuming for Argon 354 kJ/kg) 


mixture 




Ideal Argon liquefaction power 


817kW (assuming for Argon 478 kJ/kg) 


Assumed efficiency 


5% 


Estimated total plant power 


30 MW 



7.4.2 Operation 

Even with a passive insulation, the favorable area to volume ratio of the con- 
sidered tanker geometry limits the evaporation of the liquid argon to 0.03% 
of the volume per day. Hence, it would take 9 years to evaporate completely 
the tanker by evaporation (see Table EI). Once the tanker is cold and ffiled 

■^However, given the large quantities produced there is probably no market for it! 
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with argon, it will remain therefore in very stable conditions. In order not 
to loose mass, one must essentially "refill" the amount of liquid that is lost 
by evaporation. This corresponds to a loss of 0.3 liters per second. These 
requirements are modest compared to those of the initial filling!. The problem 
lies therefore in filling the tanker in a short time and not to keep it cold. In 
order to produce new liquid argon (assuming the existing one is simply venti- 
lated away, a conservative assumption), one reaches with an efficiency of 5% 
a power of 6.2 MW. If the liquid Argon is circulating in closed loop and we 
assume that the liquefaction power is still fully needed (conservative), then 
the power is 3.6 MW. These figures are summarized in Table [7| 



Table 7: Cryogenic parameters: refilling (refrigeration) 



Liquid Argon refilling rate 


0.3 liters/second or 23000 liters/day 


Argon gas equivalent 


0.2 m^/s 


Air volume equivalent (Ar 1%) 


20m Vs or 20000 1/s 


Ideal power of separation of argon mix- 
ture 


130kW (assuming for Argon 354 kJ/kg) 


Ideal Argon liquefaction power 


180kW (assuming for Argon 478 kJ/kg) 


Assumed efficiency 


5% 


Estimated total power 


6.2 MW 



7.5 Underground operation — the cavern 

We believe that such a detector should be located underground in order to 
provide the best possible physics output given its mass. With a shallow depth, 
cosmic muons induced reactions which will provide dangerous backgrounds 
(mostly neutrons from the rock) to rare events. We are currently studying 
two possible configurations: (1) a cavern in a mine (with a vertical access 
through a shaft) (2) a cavern in a mountain (with a horizontal access through 
a tunnel). Work is in progress and preliminary results from these studies 
should be available soon. 

7.6 Physics program 

The physics program of a 100 kton target with the capability to study rare 
events with the quality of a bubble-chamber are vast. The physics potentiali- 
ties of a giant liquid argon TPC with scintillation and Cerenkov light readout 
will be reported elsewhere |27j. This physics program competes favorably with 
a 1 Megaton water Cerenkov|20^. In terms of rates, one can mention: 
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• Nucleon stability: a target of 100 kton — 6 x 10^ yields a sensitivity 
without backgrounds of Tp/Br > lO^^years x T{yr) x e at the 90%C.L. 
This means that lifetimes in the range of lO^^years can be reached within 
10 years of operation. Channels like p vK have been shown to be 
essentially background free. 

• Atmospheric neutrinos: about 10000 atmospheric neutrinos per year and 
about 100 Vr charged current event per year from oscillations. 

• Solar neutrinos: about 324000 events per year with electron recoil energy 
above 5 MeV. 

• Supernova neutrinos: a SN-II explosion at 10 kpc yields about 20000 
events. 

In addition, one obtains 460 CC per 10^^ protons at 2.2 GeV (real focusing) 
at a distance of 130 km, and 15000 CC per 10^^ ^^A^e decays with 7 = 75 
at a distance of 130 km. 

8 Conclusions 

Given the tremendous physics potential of such detectors, we invite the com- 
munity to a deep reflection concerning the feasibility of giant neutrino detec- 
tors and fully compare these two options: 

• Giant 1 megaton H2O Ccrcnkov detector 

• Giant next-generation 100 kton liquid Argon detector, taking advan- 
tages of possible advances in the LAr TPC technology like a bi-phase 
operation with charge amplification for long drift distances, an Imag- 
ing-l-Scintillation-l-Cerenkov readout for improved physics performance, 
and a Giant boiling cryostat (LNG technology) 

These detectors offer the widest physics output (accelerator & non-accelerator). 
Coupled to the proper superbeams and beta-beams they could greatly improve 
our understanding of the CP-phase in the lepton sector. International sites 
with proper depths and infrastructure for potentially locating such giant de- 
tectors should be reviewed and compared. To build such large/giant detectors 
for only CP seems unconceivable, hence, giant detectors must have "broad" 
physics programs. Detectors should be underground (depth to be optimized 
vs backgrounds). 
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